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We present results on the Omega baryon {O,^) electromagnetic form factors using Nf = 2 + 1 
domain-wall fermion configurations for three pion masses in the range of about 350 to 300 MeV. We 
compare results obtained using domain wall fermions with those of a mixed-action (hybrid) approach, 
which combines domain wall valence quarks on staggered sea quarks, for a pion mass of about 350 
MeV. We pay particular attention in the evaluation of the subdominant electric quadrupole form 
factor to sufficient accuracy to exclude a zero value, by constructing a sequential source that isolates 
it from the dominant form factors. The Q,~ magnetic moment, /^q- , the electric charge and magnetic 
radius, {'"1;o/a/i)> ^^e extracted for these pion masses. The electric quadrupole moment is determined 
for the first time using dynamical quarks. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70. Dj 



I. INTRODUCTION 



The structure of hadrons, such as size, shape and charge distribution can be probed by their electromagnetic form 
factors. The baryon, consisting of three valence strange quarks is, significantly more stable than other members 
of the baryon decuplet, such as the A, with a life-time of the order of 10~^° s. This fact makes the calculation 
of its electromagnetic form factors particularly interesting since they are accessible to experimental measurements 
with smaller theoretical uncertainties. Its magnetic dipole moment is measured to very good accuracy, unlike those 
of the other decuplet baryons. A value of /iji- — —2.02(5) is given in the PDG [1] in units of nuclear magnetons 
(/ijv). Within lattice QCD one can directly compute hadron form factors starting from the fundamental theory of the 
strong interactions. Furthermore, higher order multipole moments, not detectable by current experimental setups, 
are accessible to lattice methods and can reveal important information on the structure of the hadron. An example 
is the electric quadrupole moment, which detects deformation of a hadron state. 

In this work we calculate, for the first time, the electromagnetic form factors of the 11^ baryon using dynamical 
domain-wall fermion configurations. For the calculation we use the fixed-sink approach, which enables the calculation 
of the form factors for all values and directions of the momentum transfer q concurrently. The main advantage of this 
approach is that it allows an increased statistical precision, while at the same time it provides the full dependence, 
where = —q^. In order to obtain accurate results on the form factors we construct optimized sources for the 
sequential inversion. This is particularly important for the subdominant electric quadrupole form factor, for which we 
construct an appropriate source that isolates it from the two dominant form factors [2] • This requires extra sequential 
inversions but it is essential in order to determine the electric quadrupole form factor to good accuracy. 

The form factors are calculated using Nf = 2 -f 1 dynamical domain-wall fermion configurations at the three 
lowest pion masses currently available, namely ttZtt = 350 MeV, m^r = 330 MeV and — 297 MeV. The results 
are compared to those obtained with a hybrid action that uses domain wall valence quarks on staggered sea quarks 
simulated by the MILC collaboration 

The paper is organized as follows: In Section II we provide the definitions of the corresponding multipole form 
factors and describe the lattice setup to extract them. In Section III we discuss the results and in Section IV we give 
the conclusions. 
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II. LATTICE TECHNIQUES 
A. Electromagnetic matrix element 



The n~ has spin and isospin 3/2 and therefore the decomposition of the electromagnetic matrix element is the 
same as that of the A. The on-shell fi^ matrix element of the electromagnetic current V^, is decomposed in terms 
of four independent Lorentz covariant vertex functions, ai{q^), 02(9^), ci(g^) and 02(9^), which depend only on the 
squared momentum transfer = —Q^ = {pi — PfY' ■ The initial and final four-momentum are given by pi and pf, 
respectively. In Minkowski spacetime these covariant vertex functions are given by [3] 



{n{pf,sf)\v^^\n{p,,s,)) = 



En{Pf)En{Pi) 



a2{q^) 
2mo 



q^ 
4ml 



ci{q^)r + 



C2{q^) 
2mo 



(1) 



(2) 



The rest mass and the energy of the particle are denoted by and respectively. The initial and final spin- 
projections are given by Si and s/, respectively. Recall also that every vector component of the spin-| Rarita- 
Schwinger vector-spinor u„ satisfies the Dirac equation, {pfi'y^ — mQ)u"'{p, s) = 0, along with the auxiliary conditions: 
ja-u'^ip, s) = and PcrU'^{p, s) = 0. Additionally, the covariant vertex functions are linearly related to the (dimension- 
less) electric Geo{q^)i GE2{q^) and magnetic Gj\/i(q^), Gj\/3(q^) multipole form factors [HHj- Namely, the expressions 
relating the multipole form factors and the covariant vertex functions are given in Ref. |4] and are quoted below for 
reference: 



where the positive quantity t = — 
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(5) 

(6) 
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B. Lattice setup 



We use gauge configurations generated by the RBC-UKQCD collaborations using Nf = 2 + 1 domain-wall 
fermions [F and the Iwasaki gauge-action. The simulations are carried out on two lattices of size 24^ x 64 at a 
pion mass of 330 MeV and 32'^ x 64 at pion masses of 355 MeV and 297 MeV, respectively. The latter has a smaller 
lattice spacing and therefore we will refer to it as the fine lattice. For the 24"^ x 64 lattice, or coarse lattice, the lattice 
spacing a, the light u- and d-quark mass as well as the strange quark mass were fixed by an iterative procedure using the 
n~, the pion and the kaon masses [5] as inputs. The value obtained for the lattice spacing is — 1.729(28) GeV [S]. 
For the fine lattice the scale was fixed from the ratio of the pion decay constant, calculated on the fine lattice to the 
one computed on the 24'^ x 64 at the same values of the ratio m^//^. The value found is = 2.34(3) GeV [S]. In 
addition, to these two lattices, we perform the calculation using a mixed-action with domain-wall valence quarks and 
staggered sea quarks. The gauge configurations were produced by the MILC collaboration [3 [5] using two degenerate 
fiavors of light staggered sea quarks and a strange staggered sea quark fixed to its physical mass. The lattice size is 
28"^ X 64 and the mass of the light quarks corresponds to a pion mass of 353 MeV. The lattice spacing is 0.124 fm as 
determined from the T' — T mass difference ffl. For the valence quarks we use domain wall fermions (DWF). The 
valence strange-quark mass was set using the Np — 3 ensemble by requiring the valence pseudoscalar mass to be 
equal to the mass of the Goldstone boson constructed using staggered quarks [S] . Similarly the light quark valence 
mass is tuned by adjusting the DWF pion mass to the taste-5 staggered Golstone boson pion. The domain wall quark 
masses take the values given in Table |l] Technical details of this tuning procedure are given in Refs. [9HTT] . 

In all cases we used = 16, which is what was used in the simulation of the dynamical domain wall fermions. We 
note that for the coarse lattice at the pion mass used here the residual mass is large compared to the bare quark mass 
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and chiral symmetry breaking is expected. The value of = 16 is also used in the mixed action calculation where 
it was shown that the residual mass is 10% of the bare quark mass, ensuring small chiral symmetry breaking jlO] . In 
Table H] we provide details of the simulations, along with the value of the mass of the fl^ obtained in this work as well 
as the value computed by other groups when available. 



Li 


X Lt 


' conis 


A^c^onfe- a-' [GeV] m„,d/m. [GeV] 


mjv [GeV] 


mn [GeV] 


Nf = 2 + 1 domain wall fermions [Sj 


24-^ 


X 64 


200 


200 1.729(28) 0.005/0.04 0.329(1) 


1.154(7) T2 


1.77(3) (1.758(9)) [5] 


Np = 2 + 1 domain wall fermions 6 


32'=' 
32^ 


X 64 
X 64 


200 


105 2.34(3) 0.006/0.03 0.355(6) 
120 2.34(3) 0.004/0.03 0.297(5) 


1.172(21) 
1.109(21) 


1.79(4) 
1.76(2) 








Mixed action 1181 










DWF valence: am^^a = 0.0138, am^ = 0.081 






28^ 


X 64 


210 


120 1.58(3) 0.01/0.05 0.353(2) 


1.191(19) 


1.78(3) (1.775(5)) [9| 



TABLE I; Parameters used in the calculation of the form factors. We give the number of configurations A'^^onfs used to extract 
the subdominant electric quadrupole form factor Ge2, as well as the number of configurations used A'^^o°^g to extract the 
dominant form factors for the various lattices employed in this study. The Q~ hyperon mass as determined in this work is 
given in the last column and it is compared with the value determined by the RBC-UKQCD collaboration and the LHPC for 
the mixed action as given in parenthesis. 



C. Interpolating fields 

In order to calculate the on-shell matrix element we utilize appropriate two- and three-point correlation functions. 
An interpolating field operator with the quantum numbers of the baryon is given by 

X.o(a:)=e'^''^sS(s-^''[C7.]^^s;), (8) 

where C — 7472 is the charge-conjugation matrix and a represents the vector index of the spin-| spinor. To ensure 
ground state dominance at the shortest possible Euclidean time separation we perform a gauge invariant Gaussian 
smearing on the strange quark fields that enter in the interpolating field, as described in Refs. [I3l I14j : 

s^(t,f) =^[l + ai/(f,y;;7)]""- (9) 

V 

3 

H{x, mU)^Y. -^Ul^x- fi, t)63,y-+p,) , (10) 

where s is the smeared s-quark field. The links U^ix^ t) entering the hopping matrix H are APE-smeared gauge fields. 
In particular, for DWF on the coarse lattice we have used the Gaussian smearing parameters a = 5.026 and nw = 40, 
while for the fine lattice the corresponding smearing parameters are a = 7.284 and n\Y — 84. These are the same 
parameters as those used to ensure optimal filtering of the nucleon state [6j . 

In Fig. [1] we show the results for the Qr effective mass calculated from the two-point function ratio am^^^^(t) — 

— log[G(< + 1, C))/G(i, 0)] for the three different sets of configurations considered in this study. The results are sum- 
marized in Table [U 

For the DWF simulations, on both the coarse and fine lattices considered in this work, the resulting values for the 
mass are 1.77(3) GeV and 1.76(2) GeV, respectively. These values agree with the value found in Ref. [5]. The 
same agreement is obtained in the case of the hybrid action. In Ref. [5] it was found that, at the chiral limit, the Vl~ 
mass decreases by about 2% its value at aniu,d = 0.005. Compared to the experimental value of 1.672 GeV [T] the 
value obtained at the physical point is about 50 MeV higher indicating that the strange quark mass is a few percent 
larger than the physical one in these simulations. 



4 



1.6 



1.4 -S 



1.2 



I \ 

DWF, m, = 355 MeV □ 

HYBRID, m, = 353 MeV O 

DWF, m, = 330 MeV T A 

DWF, m, = 297 MeV □ 




0.8 



0.6 



t/a 

FIG. 1: The Q~ effective mass and the fit to a constant 
plotted against the time separation for each ensemble 
considered. The statistics used to extract the effective 
masses are summarized in Table U] 
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FIG. 2: The ratio R = RaTti(r,q,t) extracted for tem- 
poral source-sink separations tf/a = 8 and tf/a = 10, 
using 50 gauge configurations. The results for tf/a = 10 
are shifted to the left by one unit. We show results for 
current direction fi — 1 and fi — 2,3 and momenta q: 
(0, 1, 0)^ and (1, 0, 0)^, respectively. The bands corre- 
spond to the constant form fit errors. 



D. Two- and three-point Correlation functions 



The electromagnetic form factors can be extracted in lattice QCD by constructing appropriate combinations of two- 
and three-point correlation functions. The corresponding lattice correlation functions are given by 

(11) 
(12) 



C(^^p,^) = 5]e-'-~/-^^r^,„ (x^c.(t,f/)xra'(o,o)) , 

Xf 

G<,^,(^^g,^) = ^e'-«0„ (x.a(t/,x/)y^(t,f)x™'(0,0) 



X,Xf 



For our lattice setup we take a frame where the final il^-state is produced at rest i.e. pf — 0. Furthermore lattice 
calculations are carried out in a Euclidean space-time, and hence from here on all expressions are given with Euclidean 
conventions |15j . We use the local vector current carrying a momentum q = —pi, which is inserted at time t. The 
renormalization constant Zy is determined by the condition Ge{0) = —1- The F matrices are given by 



r^ = ^(i + 7'), 



r'= = ^(l+7')757fe, k = 1,2,3. 



By inserting into the correlation functions a complete set of energy momentum eigenstates 



E 



VE., 



■n[p) 



\n{p,0){n{p,0\=^. 



(13) 



(14) 
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with ^ denoting all other quantum numbers, such as spin, one finds that the leading contributions for large Euclidean 
times t and tf — t are 



G,r(r^P,i) 



Mg 

En{p) 



|2-|2 g--Efnp) t i-j. [r''Af^(p)] + excited states , 

|^|2 e-^-(^'.)* tr[r^Af,,b;)Of.,,,Af,,(K)] 



excited states . 



(15) 



(16) 



The leading time dependence and unknown overlaps of the f2~ state with the initial state Jn\0) in the three-point 
correlation function can be canceled out by forming appropriate ratios that involve both the two- and three-point 
functions. The ratio employed in this work is given by the following expression 



Gkk{r\ 6,tf) V Gfcfe(r4, 6,tf - t)Gkk{T^p^,t)Gkk{r^p^,tf) 



(17) 



where a summation over the repeated indices k {k = 1,2,3) is understood. This ratio becomes time-independent 
(displays a plateau) for large Euclidean time separations, that is 



(18) 



C = 



2En{q) , 2El{q) , EUq) , EtM 



(19) 



It is understood that the trace acts in spinor-space, while the Rarita-Schwinger spin sum, expressed in Euclidean 
space, is given by 



A<Tr(p) = ^ U„{p, s)Ur{p, I 



-if + nin 
2mn 



2paPT 

3™^ 



^ Palr - Prla 



(20) 



The electromagnetic form factors are extracted by fitting R^^-^iT , q, t) in the plateau region determined by lif'^iT, q). 

Since we are evaluating the correlator of Eq. (12) using sequential inversions through the sink [16j , a separate set 
of inversions is necessary for every choice of vector and Dirac-indices. The total of 256 combinations arising from the 
vector indices of the 17" and the choice of F matrices, as can be inferred from Eq. (12), is beyond our computational 
resources, and hence we concentrate on a few carefully chosen combinations given below: 

4'^'?) = E ^m^,Av\q) (21) 

],k,i=i 

5i{En+Mn)C 



Gmi 



18M2 



[^iAl3 - 92) + S2,f,{qi - qs) + ^3,^.(92 - qi)] 



n(^\q) = Y,Uk,k{T\q) 



(22) 



fe=i 

= -Geo SM^ ^^^^^ + ^""^^''^ + "^^^^ " '^^^''^^ 
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(23) 
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QMl{En + Mn)^^ 
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^Sk,t,qiq2q3 ( 2 



gl + g2 + g3 - gfc 



M3 



(16£;o + 14Mo)gig2g3 



-10Mn(glg2 + g2g3 + g3gl)gfe 
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where the kinematical factor C is given in Eq. (19). As expected, current conservation g^IIp = is manifest in 
the right hand side of the equations. From these expressions ah the multipole form factors can be extracted. For 
instance Eq. (21 1 is proportional to Gmi, while Eq. (23) isolates Ge2 for /i = 4. Furthermore, these combinations 



are optimal in the sense that all momentum directions, each of which is statistically different, contributes to a given 
Q^-value. This symmetric construction yields a better estimator for the ri~-matrix elements than methods where 
only one momentum-vector is accessible. 

In this paper, we consider only connected contributions to the three-point function. These are calculated by 
performing sequential inversions through the sink, which necessitates fixing the quantum numbers of the initial and 
final states as well as the time separation between the source and the sink. The optimal combinations given in 
Eq. (21) - Eq. (23), from which Geo, Gmi and Ge2 are determined, can be implemented by an appropriate sink 
construction which requires only one sequential inversion for each of the three types of combinations. No optimal sink 
is considered for the octupole magnetic form factor in this work. Although it can and has been extracted, the results 
exhibit large errors and are consistent with zero. We therefore refrain from presenting this specific form factor. The 
matrix element for all the different directions of g and for all four directions /i of the current can then be computed 
yielding an over-constrained system of linear equations which can be solved for the form factors in the least squares 
sense. A singular value decomposition of the coefficient matrix is utilized to find the least squares solution. The 
statistical errors are found by a jack-knife procedure, which takes care of any possible autocorrelations between gauge 
configurations. 

As already mentioned, the three-point function of the connected part is calculated by performing sequential inver- 
sions through the sink. This requires fixing the temporal source-sink separation. In order, to determine the smallest 
time separation that is still sufficiently large to damp excited state contributions, we perform the calculation at two 
values of the sink-source separation. We use tf/a = 8 and i//a — 10 for the DWF configurations corresponding to the 
coarse lattice spacing a = 0.114 fm. We compare in Fig. [2jthe results for the plateaus IlcrT^li^, q), for a few selected 
directions of the current and for low momentum q values for these two sink-source time separations. As can be seen, 
the plateau values at tf /a = 10 are consistent with the smaller time separation the latter exhibiting about half the 
statistical error. We therefore use tf/a = 8 or tf = 0.91 fm as source sink separation. For the fine DWF lattice the 
inversions were performed for tf/a = 12, which corresponds to about tf — 1.008 fm. Similarly for the hybrid scheme 
the time separation was taken to be a,t tf/a — 8 or tf = 0.992 fm. 



III. RESULTS 



We use the local electromagnetic current, = — ^sj'^s, which requires a renormalization factor Zy to be included. 
The vector current renormalization constant is determined from the lattice calculation by the requirement that 

ZvGeo{0) = -1, (24) 



where -1 is the charge of U~ . The values of Zy extracted using Eq. ( 24 1 are given in Table |III[ where the errors shown 
are statistical. For the coarse lattice with DWF, the value of Zy = 0.7161(1) is calculated [17] from the pion decay 
constant. For the fine lattice Zy was fixed using the nucleon electric form factor [6] with values Zy = 0.7468(39) at 
= 297 MeV and Zy = 0.7479(22) at = 355 MeV. For the mixed-action HHj with = 353 MeV the value 
of the current renormalization constant Zy = 1.1169 is obtained by dividing the unrenormalized isovector current 



with the forward matrix element. These values differ by about l%-2% with the ones found using Eq. (24 1. This 
discrepancy indicates systematic errors on the 2% level. 



A. Electric charge form factor 

Our results for the electric charge form factor, Geo{Q^), are depicted in Fig. [s] for the fine and coarse lattice 
using DWF and for the mixed action. Results using the mixed action have consistently smaller values. This can be 
attributed either to cut-off effects or to a small dependence on the mass of the light sea quack mass. In order to check, 
we perform a calculation using DWF at = 355 MeV on the fine lattice for the magnetic dipole form factor. This 
will be discussed in next section. In Fig. [3] we show fits to a dipole. As can be seen, the momentum dependence of 
this form factor is adequately described in all cases by a one-parameter dipole form 



1 



(25) 



In the non-relativistic limit the slope of the above dipole form evaluated at momentum transfer 
the electric charge mean square radius by 



6 



Gbo(O) dQ 



0, is related to 



(26) 




From the dipole fit to the coarse DWF lattice data we determine Aeo and obtain a value of (r^o) — 0.353(8) fm^, 
while for the fine DWF lattice the corresponding value turns out to be (t'^q) — 0.355(14) fm^ ^. These values are 
slightly greater in magnitude than the one reported in Ref. [IS], which was obtained in a quenched lattice QCD 
calculation. The discrepancy may originate from unquenching effects or pronounced light quark mass dependence 
since the pion mass used in the quenched study of Ref. [19 is larger than what used here. The results for the (t^q) 
are given in Table |lTT] 



Note the different sign as eompared to Ref. I19| since we liere divide by Ge{0) = 



-1. 
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B. Magnetic dipole form factor 



In order to check for cut-off effects we perform a comparison between the hybrid results and results obtained at 
the same pion mass using DWF on our fine lattice. This comparison is shown in Fig. |4] The results using a hybrid 
action show a small deviation having a smaller slope as compared to the DWF results. This is the same behavior as 
was observed in the case of the electric form factor. Given the fact that the lattice spacing for the mixed action is the 
largest this points to cut-off effects. In Fig. [5] we show results obtained using DWF on the coarse and fine lattices, 
which are in agreement. This indicates that for these lattice spacings cut-off effects are small. 
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FIG. 4: The magnetic form factor Gmi{Q'^) comparing 
tlie results from the mixed action approach and the DWF 
lattice at ~ 350 MeV. 
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FIG. 5: The magnetic dipole form factor, Gmi, using 
DWF at = 353 MeV, = 330 MeV and = 
297 MeV. These results are shown along with the dipole 
fit as given in Eq. ( 28 1 . The datum for the magnetic 
dipole form factor at = GeV^ G^^p (O) = -3.60(8), 
is also included. 
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The Q^-dependence of the form factors, as in the case of G^Oj can be described by a dipole form as can be seen in 
Figs. |4] and [Sj Fitting to the two-parameter exponential, dipole and tripole forms 

GMi(g') = Goexpf-^), (27) 



GMim = (28) 

Gmi(Q')= ^ ^° (29) 

we can obtain a value for the anomalous magnetic moment of the $7". 

By utilizing the fit parameter, Gq = Gj\/i(0), and the lattice computed $7" mass from Table [ll{ we can evaluate 
the magnetic moment in nuclear magnetons, via the relation 

^in- = Go (-^\ - Go f ^iN■ (30) 



Our value of /i^- in nuclear magnetons /xjv is given in Table [lTl| The values obtained are in accord with two other 
recent lattice calculations [19l [20] . The calculation in Ref. [IQ^ is similar to ours in the sense that the three-point 
correlation function is also calculated, but the evaluation is carried out in the quenched theory and only at one value 
of Q^. In Ref. [20] a background field method was employed, where energy shifts were computed using Np = 2 + 1 
clover fermions at pion mass of 366 MeV on an anisotropic lattice. 
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FIG. 6; The subdominant electric quadrupole form factor Ge2{Q^) for DWF using the fine lattice at = 297 MeV, and 
using the hybrid action at — 353 MeV. The extrapolated values at = are also depicted. The two results, apart from 
being consistent within errors indicate a non-zero deformation for the Q~ baryon. 



type of fit Abo [GeV] 


XBo/d.o.f 


Go 


Ami [GeV] xin/d.o.l 


Nf = 2 + 1 DWF (24^ X 


64), iVconfs= 200 








exponential 

dipole 1.151(13) 
tripole 


1.500 


-3.264(89) 
-3.601(109) 
-3.478(101) 


1.016(27) 
1.187(41) 
1.555(49) 


0.225 
0.860 
0.435 


Nf = 2 + 1 DWF 


(32^ X 64), iVconfs = 105; 


= 355 MeV 




exponential 
dipole 
tripole 




-3.246(96) 
-3.557(130) 
-3.443(116) 


0.996(43) 
1.171(63) 
1.530(76) 


0.159 
0.440 
0.240 


Nf = 2 + 1 DWF 


(32^ X 64), iV,o„fa = 120; 


= 297 MeV 




exponential 

dipole 1.146(23) 
tripole 


0.887 


-3.199(155) 
-3.443(173) 
-3.355(165) 


1.061(48) 
1.277(68) 
1.656(83) 


0.080 
0.064 
0.040 


Hybrid (28"* x 


64), iVconfs = 120 








exponential 

dipole 1.213(17) 
tripole 


0.168 


-3.154(69) 
-3.368(80) 
-3.293(76) 


1.064(30) 
1.285(44) 
1.662(54) 


1.147 
0.163 
0.053 



TABLE II: The fit parameters for the exponential, dipole and tripole forms extracted from the lattice data. For the fine lattice 
with = 355 MeV DWF we have only performed inversions for the source type associated with the dominant magnetic dipole 
form factor GMiiQ^) (see Eq. 



C. Electric quadrupole form factor 



From the perspective of hadron structure, the extraction of the electric quadrupole form factor is of special interest 
since it can be used to provide valuable information regarding the deformation of a hadron. In this work we extract 
for the first time in unquenched QCD the subdominant Ge2 form factor for the fl^ baryon, to sufficient accuracy to 
exclude zero values. This has been achieved by utilizing two different lattices: namely, the fine DWF lattice and the 
MILC lattice at lattice spacings of a = 0.084 fni and a — 0.124 fm respectively. We note that for the coarse DWF 
lattice the results for Ge2 are too noisy to exclude a zero value and we therefore do not present them here. The lattice 
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lattice TT^TT Zv 






Q 3 

2 






[Ll X it] [GeV] 


[mn] 


[fm^] [fm^] 




This work 


HYB: 28^ x 64 0.353 1.121(2) 


-1.775(52) 


0.283(20) 0.338(9) 0.838(19) 


-1.366(222) 






DWF: 24^ X 64 0.330 0.727(1) 


-1.904(71) 


0.332(23) 0.353(8) 


— 






DWF: 32^ X 64 0.355 0.7479(22) 


-1.868(78) 


0.341(37) 








u vv r . X 04 u.zy( V. f o^oy^) 


-1.835(94) 


0.286(31) 0.355(14) 0.959(41) 








extrapolated 0.140 - 


-1.875(399) 0.321(16) 0.348(52) 0.898(60) 


-1.651(262) 


Ref. 




20^ X 40 0.697 1 


-1.697(65) 


0.307(15) 




Ref. 




24^ X 128 0.366 


-1.93(8) 






Ref. 


1 




-2.02(5) 







TABLE III: The magnetic moment , the electric charge and magnetic radii and the electric quadrupole moment as 



extracted using Eq. |3l|. The values of 



fi)i ('"bo) ''■nd Q§ shown above arise from the dipole fit form. Note that 



rln) 



_6 dGMl(Q^) I 



l(0) 



dQ2 



results for Ge2 are depicted in Fig. [6j The value of the quadrupole electric form factor Ge2{Q^) at = using the 
exponential form to fit the lattice results is 0.756(298) for the hybrid action and 0.882(475) for the fine DWF lattice. 
From these results it is readily deduced that the shape of the fi~ hyperon must deviate from the spherical one. 
The electric quadrupole moment determined from the fits as Qci = G£;2(0)^^ can be related to the transverse 

charge density in the infinite momentum frame. For instance, the transverse charge density defined in the light-front 
for spin projection 3/2, is given by pTl 



g? = -h[GMi(0)-3eo] + [G£2(0) + 3eo] 



(31) 





FIG. 7: Transverse charge densities in the with polarization along the x-axis. Left: Pt3/2(6)- Right:pyiy2(^)- ^ circle of 
radius 0.5 fm is drawn in order to clearly demonstrate the deformation. For the evaluation of the densities we used the dipole 
parametrization of the form factors. 



We note that for a spin-| particle without internal structure, for which Gi\ii{0) = 3en and Ge2{0) = — 3eo pTll22j . 
the quadrupole moment of the transverse charge densities vanishes. We calculate this qua ntity by using a fit to the 



electric quadrupole to obtain the value at = 0. The results obtained are shown in Table |III| and plotted in Fig. [10 
for the dipole fitting Ansatz. Both of the two values are negative and consistent within statistical errors. Therefore, 
they suggest that the quark charge distribution in the il" must be deformed. In order to investigate the deformation 
in more detail we construct the transverse charge density in the infinite momentum frame, following Refs. [2T | 122 ) . 
Considering the spin of the along the x-axis and states of transverse spin s± =3/2 and Sj^ = 1/2 we obtain the 

transverse charge densities p^^^^^^) and in term of the two-dimensional impact parameter h. In Fig. jTj we 

compare and p^^^^^)- can be seen, in a state of transverse spin projection s±^ — 3/2 the shows a 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

by I fm] 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

by Ifm] 



FIG. 8: Comparison of the transverse charge densities Pr3/2(^) (left) and PTi/2ib) (right) along the y-axis to the monopole-field 
(symmetric) shown by the dashed line. 

small elongation along the spin axis (prolate) ^. This elongation is less as compared to that seen for the A+. As in 
the case of the A+, in a state of transverse spin projection = 1/2 the is elongated along the axis perpendicular 
to the spin. 




-1.5 



FIG. 9: The individual multipoles for the transverse charge density Pr3/2(^) in the Q, with polarization along the x-axis. Upper 
left: monopole-field. Upper right: dipole-field. Lower left: quadrupole-field. Lower right: Octupole-field 



^ Note that this is consistent with the negative sign of Q3/2 since the f2 is negatively charged and have included its charge in the 
electromagnetic current 



12 



In Fig. [8] we show the profile of the transverse densities compared to the monopole field that is symmetric. In Fig. [9] 
we show the individual multipole fields for the state with transverse spin s± =3/2. 



D. Extrapolation to the physical point 



In this section we examine the sea quark dependence of the magnetic moment, radii and the quadrupole moment. 



They are extracted by fitting the Q^-dependence of the form factors to a dipole form. As can be seen from Fig. [10 
the sea quark mass dependence is consistent with a constant for all quantities confirming that sea quark effects are 
small. In particular, the value of the magnetic form factor at = is consistent with experiment. On the other 



hand extrapolating the magnetic moment we obtain the value given in Table III This is 5% smaller than experiment 
which is to be expected given the larger value of the strange quark mass. The reason is that the mass of the fi" is 5% 
larger than experiment and this will affect the value of the magnetic moment when we convert to nuclear magnetons. 
In the fits for the magnetic moment and radii we did not include the results obtained in the hybrid action because of 
the small finite a effects observed. Given the large statistical errors on quadrupole moment such small finite-a effects 
are negligible and therefore, in this case, we include the result using the hybrid action to obtain the extract the value 
at the physical point. In Table III we give the values that we find at the physical point for the radii and the dipole 
and quadrupole moments of the transverse charge density obtained from Eq. (31 1. 



a 



-3.2 
-3.4 
-3.6 
-3.8 

.4 

.3 

.2 

.1 

.0 

.38 

.36 

.34 

.32 

.30 
-1.2 

-1.4 



-2.0 



O Hybrid: m,=0.353 GeV 
A DWF: m =0.330 GeV 



n DWF: fine m,=0.355 GeV 
A DWF: fine m =0.297 GeV 




0.00 



0.05 



0.15 



0.10 
(GeV^) 

FIG. 10: From top to bott om we show Gmi{0), the magnetic radius {r\ji), the electric radius {r%Q) and the quadrupole 
moment extracted from Eq. (31 1 as a function of extracted from dipole fits. The point shown by the filled square is the 
value extracted from the fit at the physical pion mass. In all cases except for the quadrupole moment the results using the 
hybrid action are excluded form the fit. 
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IV. CONCLUSIONS 



By utilizing properly constructed sequential sources the dominant il^ electromagnetic form factors Geo and Gmi 
are calculated with good accuracy using dynamical domain-wall fermion configurations as well as a hybrid action. 

In addition, we extract the magnetic moment of the fl^ by fitting the magnetic dipole form factor Gmi to a two- 
parameter dipole form. We find a value that is within errors to the experimentally measured value [T]. The electric 
charge and magnetic radii {{r\f^) and {r\ji)) are computed and like the magnetic dipole moment they do not show 
sea quark dependence in the range of masses studied in this work. 

Finally, the subdominant electric quadrupole form factor G e2 is computed for the first time in an unquenched lattice 
calculation to sufficient accuracy to exclude a zero value. This has been accomplished by constructing an appropriate 
sink that isolates it from the two dominant form factors. We find consistent results with DWF and using a hybrid 
action. The positive non-zero values of Ge2 Q'^ — ^ suggest that the structure of the baryon is non-spherical. 
In the light-front frame we find that the quark charge density in a VL~ for a state of transverse spin projection -1-3/2 
is shows an elongation along the axis of the spin (prolate deformation). As compared to the A"*" in the same state 
the amount of deformation seen in the 57" is smaller. 
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